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Recent theoretical developments to calculate cross sections of liadronic objects in the high energy limit 
are summarised and experimental attempts to establish the need for new QCD effects connected with 
a resummation of small hadron momentum fractions x are reviewed. The relation between small-x 
parton dynamics and the phenomenon of diffraction is briefly out-lined. In addition, a search for a 
novel, non-perturbative QCD effect, the production of QCD instanton induced events, is presented. 



1 Mini-Introduction to Recent 
BFKL Progress 

When two hadronic objects collide at high 
energies, /ii /12 h^hiX, often more than 
one energy scale is involved. Of particular 
interest are collisions where the squared cen- 
tre of mass (cm) energy s = {hi -t- /i2)^ is 
much larger than the squared transverse mo- 
mentum transfer t — (hi — h^)^ , i.e. s ^ t ^ 
^QCD- In perturbative QCD (pQCD) calcu- 
lations of such processes large logarithms of 
the ratio of these scales, (a^ In (s/t))", arise 
at each order n in the strong coupling con- 
stant as. It is usually assumed that they 
must be resummed to obtain a reliable cross 
section prediction. 

About 25 years ago, Balitsky, Fadin, Ku- 
raev and Lipatov (BFKL) derived in the 
leading logarithmic approximation (LLA) an 
equation describing the resummation of a cer- 
tain class of exchanges with orte chain of mul- 
tiple gluons in the i— channefl. It predicts a 
power-law rise of the scattering cross section 
at high energies: 
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where is taken at the value of a hard 
scale involved in the process. Since for typ- 
ical values of as, e.g. as = 0.2, A ~ 0.5, 
the predicted rise of the cross section is so 
strong that the unitarity bound will be vio- 
lated. This is an indication that at high en- 
ergies novel effects like parton recombination 
or multiple perturbative scatterings taming 
down the rise should be important. Over the 



Figure 1. Sketch of a high energy collision of two 
hadronic objects hi /12 — > h^, X. 

past years much experimental effort has been 
devoted to observe to such phenomena. Al- 
though some hints for a rise of the scatter- 
ing cross section at high energies have been 
found, the predicted large value of A could 
not be confirmed. 

After a huge theoretical effort by many 
authors, the calculations of the next-to- 
leading order (NLO) corrections ±0 the BFKL 
equations were finalised in 19980. The NLO 
terms are suppressed by a power of 
as{as In (s/t))". To the consternation of the 
community, the NLO corrections turned out 
to be large and lead to much smaller values of 
A than the leading order (LO) contribution. 
Fig. 1^ shows A calculated in LO and NLO as a 
function of as. For instance, for a scale where 
as ~ 0.16 the NLO corrections exactly can- 
cel the LO term. For larger values, A even 
becomes negative. Naively, this would lead 
to a cross section that decreases rather than 
increases as a power of s. This is in contra- 
diction to the data. The large NLO correc- 
tions put the applicability of the BFKL re- 
summation in question. It seems necessary to 
understand the physical reason for the large 
corrections and to include them at all orders. 
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Figure 2. The BFKL exponent A as a function of 
as in the LO and NLO approximation and the 
coUinearly enhanced resummation included at NLO. 

Several solutions have been proposed to 
solve this problemaa, all of which are related 
to modelling and to resumming the higher 
order contributions. One promising method 
is to analyse the structure of the divergences 
of the BFKL characteristic function^ at all 
orders by studying thejimit where partons 
are emitted collinearlyQ. The collinearly- 
enhanced contributions have been shown to 
be able to estimate the NLO corrections and 
there are good reasons to believe that they 
also might give a significant contribution of 
the higher order corrections beyond NLO. 
The exponents obtained from this coUinear 
resummation are shown in Fig. ^. The label 
LOs is related to the power of the cross section 
expected for 7*7* collisions or jet observables 
in collisions of hadronic objects. The line la- 
belled LUc is the leading pole in the anoma- 
lous dimension of the gluon. The difference 
between the two exponents does not reflect 
an additional uncertainty, but is caused by 
additional corrections for quantities involv- 
ing an effective cut-off on the lowest accessi- 
ble momentum necessary to avoid diffusion 
into the infra-red and non-perturbative re- 
gion and is closely related to differences in 
running coupling effects. The resummation 



" The characteristic function is the Mellin transform 
of the BFKL kernel. The anomalous dimension is th« 
Mellin transform of the splitting function. See ref|3 
for an introduction to more details. 




Figure 3. Sketch of small-x signatures in high energy 
collisions: a) total cross section in DIS, b) dijet pro- 
duction in pp collisions, c) total hadronic cross section 
in 7*7* collisions d) forward jet production in DIS. 

in the coUinear limits leads to values of A 
which are in qualitative agreement with the 
experimental observations (see later). How- 
ever, a quantitative analysis is still lacking. 

Although most of the theoretical effort 
was devoted to understanding the BFKL be- 
haviour at NLO, most of the phenomenolog- 
ical studies still use LO. For instance, only 
recently part of the virtual photon impact 
factor, i.e. the coupling of the photon toJJie 
gluon chain, has been calculated in NLOlla. 

2 Parton Evolution and DIS 

In deep-inelastic scattering (DIS)[|, a pho- 
ton with virtuality S> rrip, where nip 
is the proton mass, collides with a proton. 
When the squared j*p cm energy is much 
larger than Q^, the momentum fraction of the 
struck quark, x « Q'^/W'^, becomes small. 



''For an introduction and literature on DIS, DGLAP, 
small-x and related topi cs see: 



|ittp:/ /sesam.desy.de/'-'^ carli/hcra.html 



In this case, the simple picture of DIS as a 
process where a virtual photon interacts in- 
stantaneously with a point-like parton quark 
freely moving in the proton has to be modi- 
fied. Indeed, the probability that a gluon ra- 
diates becomes increasingly high at small- a; 
and therefore the quark struck by the pho- 
ton originates from a cascade initiated by 
a parton with high longitudinal momentum. 
A sketch of such an interaction is shown in 
Fig. ^). The gluon is the driving force be- 
hind the cascade indicated by the circle. 

There are three equations describing the 
parton evolution in different kinematic re- 
gions: DGLAP, BFKL and CCFM equations. 

2.1 DGLAP and DLL Evolution 

To correctly calculate the inclusive DIS cross 
section (for all x- values), terms of the form 
{as In (Q^/Qo))" have to be resummed, since 
the smallness of as is compensated by the 
large size of In Q^. This resummatiojQ can be 
achieved by the DGLAP equationsQ. They 
are related to the renormalisation group 
equation and describe the change of the par- 
ton density in the proton with varying spatial 
resolutions . Once the parton distribution 
is specified at a given scale Qq, it can be pre- 
dicted for any other . In an appropriate 
gauge, the subsequent parton emissions from 
the proton to the struck quark (connecting 
the soft proton constituents to the hard sub- 
process) can be interpreted as a ladder dia- 
gram whose rungs and emitted partons are 
strongly ordered| in transverse momentum 
kr- This leads to a suppression of the avail- 
able phase space for gluon radiation. 

For high energy scattering, i.e. at small- 
x, the ladder becomes long and consists dom- 
inantly of gluonsj^. The correct description 

'^kT ordering is an assumption necessary to deri\f« the 
DGLAP equations. For an introduction see ref.Q. 
"^The splitting function Pgg(z) (describing the radia- 
tion of gluons from gluons) dominates at small lon- 
gitudinal momentum fractions z (Pgg(z) ~ I/2: while 
Pqq{z) ~ constant). 2 ranges between x and 1. 



of the small- a; region is still a matter of de- 
bate and subject to many technical difficul- 
ties. The dominant terms are of the form 
(a^ In (1/x) In (Q^/Qo))"- Their resumma- 
tion is referred to as the double-leading loga- 
rithmic approximation (DLL). Both longitu- 
dinal and transverse momenta are strongly 
ordered. The DLL leads to a gluon den- 
sity strongly rising with an effective power 
A w J{12as/TT) ln(l/2;) ln(Q2/<3o) towards 
low-xB. A depends on x and Q^, e.g. it in- 
creases towards higher Q^. 

2.2 BFKL Evolution 

For Ri Qq, double logarithms no longer 
dominate and the terms {as ln(l/x))" « 
{as In (W^^/Q^))" are important and must be 
resummed using the BFKL equation. In a 
physical gauge, these terms correspond to an 
rt-rung ladder diagram in which gluon emis- 
sions are ordered in longitudinal momentum, 
z. The strong kx ordering is replaced by a 
diffusion pattern as one proceeds along the 
gluon chain. The BFKL equations describe 
how a particular high momentum parton in 
the proton is dressed by a cloud of gluons at 
small-x localised in a fixed transverse spatial 
region of the proton given by Q^. 

A framework to describe the inclusive 
DIS cross section at small-x is a kx depen- 
dent, unintegrated, gluon density /{z^kx). 

/dz 
— d kx cr{x/ z, kx) f{z, kx) 

where a{x/ z^kx) is the partonic cross sec- 
tion. Together with the kx factorisationtj 
this forms a basis for a general calculation 
scheme for small- a; processes. 

At very small- a;, it is expected that many 
gluons coexist and will no longer act as free 
partons, but interact with each other. This 
"saturation" regime is characterised by an 
equilibrium of gluon emission and absorption. 
In the regime accessible for HERA such ef- 
fects are expected to be smalltS, but could 
be rather large when lower x-values could be 
reached. At THERA colliding the 500 GeV 



electrons from TESLA with the 920 GeV pro- 
ton from HERA such effects will be visible. 
Also if HERA collided heavy ions instead of 
protons, an experimental study of saturation 
effects would probably be possible. 

2.3 CCFM Evolution 



The CCFMiJ evolution equation is based on 
angular ordering and colour coherence. As 
a result in the appropriate limit it repro- 
duces the DGLAP and the BFKL approxi- 
mation. The angular ordering of the emit- 
ted partons in the initial cascade results from 
colour coherence and acts as unifying prin- 
ciple^. At small values of the parton mo- 
mentum fractions z, a random fc^ walk is 
obtained. A first attemptiJ to incorporate 
the CCFM equations in a Monte Carlo sim- 
ulation program was not suited for practi- 
cal purposes. Later a further analysisEj to 
make use of the CCFM equations for small- 
X phenomenology showed the importance of 
the inclusion of soft corrections and revealed 
many theoretical difficulties, in particular 
with terms diverging for high rather than low 
momentum fraction z. Recently, significant 
progress has been made to correctly imple- 
ment the CCFM evolution equations and to 
efficiently generate unweighted events using 
the backward evolution in Mnntp Carlo sim- 
ulation program CASCADEllaO. A slightly 
different approach is followed in the LDC 
programlla. Both programs are able to make 
predictions for any physical observable based 
on the hadronic final state. Only gluon- 
initiated processes can be simulated up to 
now. Some problems still remain, e.g. so far 
as in the CCFM evolution only the singular 
terms in the splitting functions (for z 
and z 1) have been considered. 



"^From the proton p to the photon in the coUinear 
DGLAP approximation fcy and therefore the angle 
increases, while in the BFKL approximation the an- 
gle increases because of the decrea^iig longitudinal 
momenta {9 ^ kx/izp)). See ref.E3 for an intro- 
duction to angular ordering as interpolation between 
DGLAP and BFKL. 



3 Small-x Signatures in High 
Energy Collisions 

Several signatures for small- a; QCD effects 
have been recently investigated: jets at large 
rapidity in pp collisions at TEVATRON, the 
total hadronic cross section in 7*7* collisions 
at LEP and the rise of F2 towards small- a; and 
the production of forward jets and particles 
in ep collisions at HERA. 

3.1 Jet Production at Large Rapidity in 
pp Collisions 

In the high energy limit of pp collisions inclu- 
sive dijet production (see Fig. ||b) is sensitive 
to large logarithms of the form: 



cr DC exp (A In 



) oc exp (A At]) (2) 



where ~ Et,iEt.2 is the momentum 
transfer to the hard scattering and W"^ = 
X1X2S is the squared cm energy of the hard 
scattering proces^. The last proportionality 
in equation || follows after some calculations 
for A77 ^ 1 from 

Xi = ^^li exp [rf) cosh (3) 

2 

with i = 1, 2 and 77 — (771 + rji)/2 and At; = 
I '71 ^ for the most backward (771) and the 
most forward (772) jet. Equation |^ is only 
valid for 2-body kinematics. 

The largest logarithms are therefore en- 
hanced for minimal jet Et and maximal 
W . To avoid the strong dependence on the 
steeply falling parton densities at small-.T it 
has been proposed by Mueller and Navelettj 
to study the dijet cross section at fixed Xi 
at different cm ene rgies s. In an analysis 
performed by DOcJ, jets have been selected 
by ET,t > 20 GeV and \r],\ < 3. Moreover, 
400 < (Q2 ^ Et,iEt^2) < 1000 GeV^ was 



^ Here Ex^i are the transverse jet energies, rn the 
pseudo-rapidities and Xi the momentum fractions of 
the two incoming partons. 



demanded. The cross section ratio at fixed 
Xi for two cm energies lias been measured: 

_ exp (A (Ajy^ - A??^)) 



R 



(4) 



Tlie result of the cross section ratio at = 
1.8 TeV and y/sb — 630 GeV is shown as 
a function of the mean (A?]) for y/sb = 
630 GeV in Fig. |. At large A77 > the di- 
jet cross section increases almost by a fac- 
tor of 3 between the two cm energies. The 
strong increase is reflected in a large mean 
A value: A — 0.65 ± 0.07, stronger than ex- 
pected by theoretical predictionsEJ. The ex- 
act LO pQCD calculation leads to a falling 
cross section. The LO BFKL calculation (la- 
belled LLA in Fig. |) predicts A = 0.45 for 
a^(20 GeV) = 0.17. A complete NLO BFKL 
calculation is not yet available. Surprisingly, 
the highest A = 0.6 is obtained by HERWIG, 
a LO QCD event generator supplemented by 
leading log parton showers. 

However, it has recently been pointed 
ouiE3 that an interpretation of these results 
is difficult because of differences in the defi- 
nition of the cross sections between the DO 
data and the original Mueller-Navelet pro- 
posal. This concerns the reconstruction of 
the momentum fractions Xi using eq. |3| which 
is only valid for 2-body kinematics and the 
presence of an upper bound on the momen- 
tum transfer . These differences to the 
original proposal can be neglected only in 
the asymptotic limit at large s and large Arj. 
At TEVATRON, however, this limit is not 
reached. Moreover, the requirement of two 
jets with the same minimum Et is partic- 
ularly criticalii, since large logarithms not 
connected with BFKL effects make an inter- 
pretation very difficult. 

3.2 Inclusive Hadronic Cross Section in 
7*7* Collisions 

By studying the total hadronic cross section 
7*7* — > hadrons in e~^e~ collisions, difficul- 
ties connected with the hadronic nature of 
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Figure 4. Ratio of the DO dijet cross section at a cm 
energy ^ = 1800 GeV and = 630 GeV for Ar) > 
1 and Ari > 2 as a function of the mean rapidity 
difference of the jets at = 630 GeV. 

incoming particles can be naturally avoided. 
Virtual photons are colourless objects and 
their virtuality controls the transverse 
size oc of the hard processes. For 

^QCD ^ complete perturbative calcu- 
lation is possible. For small virtualities Qf 
of one of the virtual photons and for large 
cm energies W of the 7*7* system the cross 
section contains large logarithms of the form: 



a~f*-y* oc exp (A In 



) = exp(A Y). 



If ~ Ql, a DGLAP parton evolu- 

tion between the two photons is suppressed 
and a resummation of the terms Y = 

is there- 



in (l^VVQi^?!) is needed, cr. 
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LO (Bartels 99), Q = IS GeV 
NLO(Kim99), = 15 GeV^ 



Y = ln(wVlQ,Q,l) 

Figure 5. Total hadronic cross section for the reaction 
7*7* — > X as a function of Y calculated in a LO and 
an approximate NLO BFKL calculation. 

fore often considered as a golden BFKL sig- 
nature. An example of a Feynman diagram 
is shown in Fig. ||c. As can be seen in 
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Figure 6. Total hadronic cross section 7*7* — > Jf as 
a function of the variable Y a) measured by L3 and 
b) by OPAL. Overlayed is a LO and NLO fixed order 
calculation. 



Fig.p, in LO a strnng increase of ti-y.-y* at 
high Y is expected^. An approximation to 
a NLO calculations leads to a much sup- 
pressed cross section at high Y. 



In a first measurement of a. 



7*7* 



collisions at y/s^e = 91 — 183 GeV the pre- 
dicted significant rise towards large Y was ob- 
served by L3E3. However these data have not 
been corrected for QED radiative effects. In 
particular, if the kinematics is reconstructed 
using the scattered electrons, these correc- 
tions are substantial. In a recent analysi£], 
L3 has presented cross sections corrected for 

radiative effects using e+e~ data at a cm en- 3,3 Inclusive DIS: F2 



at low Y (Y < 5). However, at the largest 
Y {5 < Y < 7), the calculations are below 
the data by 4 standard deviations. The NLO 
corrections are relatively small and they only 
slightly increase the cross section at large Y. 
L3 concludes that "this can indicate the pres- 
ence of QCD resolved processes or the onset 
of BFKL dynamics"E3. 

However, this effect has not been con- 
firmed by OPAlS analysing e~^e colli- 
sions at cm energies of ^See = 189 — 209 
GeV. The measured cross section ct^*^* (for 
Ei > OAEb, where Eb is the beam energy, 
34 < e, < 55 mrad and W > 5 GeV) 
for an average 

Q2 = 17.9 GeV^ is shown 
in Fig. Ijb) as a function of Y. The exact 
LO and the NLO calculation reproduce the 
data well. Also an approximate higher or- 
der BFKL calculationE2l imposing energy mo- 
mentum conservation along the gluon ladder, 
is in agreement with the data (not shown). 
Although in the L3 analysis slightly higher 
W and Y values are reached, the different 
conclusions which can be drawn from the 
two analyses is surprising. Some experimen- 
tal questions in particular the control of the 
radiative corrections as a potentially large 
source of systematic uncertainty needs to be 
rediscussed^ before firm conclusions can be 
drawn. To unambiguously establish small- a; 
effects in 7*7* collisions one probably needs 
much higher cm energies, like the one that 
wiU be available at TESLA. 



ergy of y^s^ = 189 - 209 GeV. The QED 
corrections lower the measured uncorrected 
cross section by 58% at large Y and, more- 
over, can be very different for the various 
contributing QCD subprocesses. The cross 
section measured in the kinematic region: 
4 < Qf < 44 GeV^ W > 5 GeV, electron 
energy Ei > 40 GeV and the polar elec- 
tron angle 30 < 0^ < 66 mrad is shown in 
Fig. ^). A fix£d order calculation in LO as 
well as in NL0c3 is able to describe the data 



A classic key measurement to understand 
the structure of the proton and its dynam- 
ical processes is the precise determination 
of the structure function F2 by counting in- 
clusively the lepton scattered ofFthe pro- 
ton. The HERA measurementsE2l'E3 span 
about 6 orders of magnitude in Q^, 0.1 ^ 



®One question is e.g. to which extent the size of the 
QED corrections depend on the Born cross section 
itself. An iterative procedure might be required. 
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Figure 7. Slope of the proton structure function A as 
a function of x in bins of . 

g2 < 30000 GeV^ and cover the momen- 
tum range from the sea quark region at low 
X (a; > 10^^) to the valence quark region at 
large x (x< 0.65). 

Compared to observables based on the 
hadronic final state, the measurement of 
has the advantage that it can be directly 
compared to pQCD calculations. However, 
because of its inclusive nature, small novel 
effects on top of the dominating DGLAP 
parton evolution might be difficult to re- 
veal. Therefore F2 has to be precisely pinned 
down to get a handle on the QCD evolu- 
tion and to simulanteously constrain the non- 
perturbative parton density functions. In 
the most recent F2 measurements a preci- 
sion of 2 — 4% in the high statistics region 
Q2 < 100 GeV^ has been achievedS. One of 
the most important results is that the rise of 
F2 towards small-x: 

1/1/ .I- 

expected in the DLL or BFKL limit has been 
observed. However, analyses based on the 
NLO DGLAP equations give a satisfactory 
description of the d a t a dow n to remarkably 
low and x valucsES'EJca. In these analy- 
ses the quark and gluon densities are param- 
eterised at a starting scale and are adjusted 
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Figures. Slope of the incl. total 7'*' p cross section A 
as a function of . Shown is a generalised vector 
dominance model and a QCD NLO DGLAP fit. 

such that the data are described. While F2 
constrains mainly the quark densities, the 
scaling violation dF2/d\ogQ^ give access to 
the gluon density. Whether this proves the 
correctness of the DGLAP picture or is just 
due to a large fiexibility of the fit, in par- 
ticular in the parameterisation of the parton 
density function, can not yet be answered. 

To improve the sensitivity ±0 small- a; 
QCD effects, the HI collaborationLj has mea- 
sured the slope dF2/dlnx for fixed as a 
function of x (see Fig. |^) . The slope does not 
change with x and the NLO DGLAP fit de- 
scribes the data. dF2 /d\nx can therefore be 
identified with the slope A only changing with 
Q^. The dependence of A is even more 
visible in the ZEUS measurement shown in 
Fig. H (see also sec. ^. Despite this success 
it is worth noting that in most NLO DGLAP 
fits the longitudinal structure function F^ be- 
comes negative, which is obviously unphysi- 
cal. The inclusion of Inll/a;) terms is able 
to cure these problemslO'^. If one aims for 
a global description of all available data, the 
standard NLO DGLAP fit has to be stretched 
to the edgeEj. In particular, there is a com- 
plicated interplay between the small- a: DIS 
data and the high- a; incliisiye jet cross sec- 
tions from TEVATRON00. The jet data 
prefer a high gluon at large x values (a; > 0.1) 
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Figure 9. Proton structure function data as a func- 
tion of at fixed x compared with several different 
fits within the MRST global fit framework. 

and this significantly influences the small-x 
gluon via the momentum sum rule. The DIS 
data prefer a smaller gluon at high-x values. 
Therefore, some kind of compromise has to 
be found. Within the present accuracy of the 
data, this is still possible, but is seems that 
there is not much freedom left in the NLO 
DGLAP fit. 

Within the framework of the MRST 
global QCD analysisEjH, it has been shown 
that problems in describing HERA F2 data at 
medium-x values 10~^) can be avoided, if 
the data at small- a; {x < 10"'^) are not con- 
sidered in the fit or, alternatively, if a small- 2: 
resummation obtained by solving the running 
coupling BFKL equation is includedEZlcj. 
This is demonstrated in Fig. where the 
standard NLO DGLAP fit gives a bad de- 
scription of the data at x — 1.3 • 10~^ and 
100 <Q^ < 1000 GeV^ while a fit only in- 
cluding data with a;>5 • 10^^ is much bet- 

''I thank R. Thorne for providing me with this figure. 



ter. This fit, however, does not describe the 
small-x data. When small-x resummation is 
included the data can be described over the 
whole kinematic domain. 

In conclusion, the NLO DGLAP fit de- 
scribes DIS data to surprisingly small- a; and 

values. However, there are some indi- 
cations that the fiexibility in the input pa- 
rameters is stretched to the edge, in partic- 
ular, if the new inclusive high Ex jet data 
from TEVATRON are included. Fits includ- 
ing small-x resummation lead to an improved 
description of the data and a better conver- 
gence of Fl when going from LO to NLO. 

3.4 Heavy Quark Production 

The production of heavy quarks Q in DIS is 
directly sensitive to the gluon content of the 
proton, since in LO the process jg — > QQ 
dominates. Heavy quark cross sections are 
therefore an excellent testing ground for the 
dynamics of the parton evolution. The in- 
clusive D* meson crosssection has recently 
been measured by H1e3 in the range 1 < 
g2 < 100 GeV^ 0.05 < y < 0.7, pf* > 

1.5 GeV and |r/^'| < 1.5 to be: ayis = 
8.5±0.42(stat)tH2(syst)±0.65(model) nb. A 
NLO calculationO based on the heavy quark 
matrix elements, DGLAP parton evolution 
and the Peterson fragmentation functiono 
only predicts ayis = 5.1 — 7 nb depending on 
the exact value of the charm quark mass rric 
and the Peterson fragmentation parameter e. 
The most significant difference between data 
and NLO is observed towards the proton di- 
rection, i.e. at large rj. Similar results have 
been obtained by ZEUS0. A possible expla- 
nation that this effect is due to neglecting the 
colour force between the charm quark and the 
proton remnantcZI is not large enough to ex- 
plain this discrepancy. The CCFM predic- 
tion, based on the off-shell matrix elements 
and on an unintegrated gluon density ob- 
tained from a £3, fit, has been calculated us- 
ing CASCADEt£l and agrees with the data. 

The visible cross section can be converted 
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Figure 10. Charm proton structure function i^'j as a 
function of logo; in bins. 

to the semi-inclusive charm structure func- 
tion F2 which is easier to compare to theory 
predictions. However, the extraction of F2 
faces the problem that only part of the total 
cross section (typically 30%) can be measured 
in the detector. This quantity is therefore 
significantly biased by the theoretical model 
used to determine the experimental accep- 
tance (up to 20% at the smallest a:). F2 ex- 
tracted in the NLO DGLAP scheme is shown 
in Fig. 10, A steep rise towards small-x and 



a clear dependence of the slope A on Q is 
observed. In particular, in the low region 
this rise tends to be steeper than predicted 
by the NLO DGLAP calculation^. At small- 
X, F2 grows at the same rate as the fully in- 
clusive structure function F2 as is expected, 
since both processes are initiated by gluons. 

This year HI has reportedca on the first 
measurement of beauty production in DIS in 
the region 2 < < 100 GeV^ 0.05 < y < 
0.7 for muon polar angles 35° < 9 < 130° 
and transverse momenta pt > 2 GeV. The 
beauty cross section is derived by determin- 
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'This is true when compared to the HI NLO DGLAP 
fit. In the ZEUS fit a slightly larger error band is 
obtained making the above conclusion less clear. 
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Figure 11. Ratio of the measured to the calculated 
total beauty production cross section in ep collisions. 



ing the fraction of beauty events decaying 
semi-leptonically to muons using the trans- 
verse momentum of the muon with respect 
to the jet axis {PT,rei) and by exploiting 
the finite lifetime of b-mesons by measur- 
ing the distance of closest approach of the 
muon track with respect to the primary event 
vertex (impact parameter). The measured 
cross section aep^bx^fiX = 39 ± 8(stat) ± 
10(syst)pb_is significantly above the NLO 
prediction^ of cr — 11 ± 2pb. In the 
calculation, a beauty mass mf, = 4.75 GeV, a 
renormalisation and factorisation scale ^ = 
\/ Q"^ + 4m^ and a Peterson fragmentation 
parameter e = 0.0033 are used. The LO 
cross section is only cr = 9 pb. The CAS- 
CADE program based on the CCFM equa- 
tion predicts cr = 15 pb. A similar ex- 
cess has been reportedin quasi-real photon- 
proton collisions by Hlc3, where the cross sec- 
tion for hb production for < 1 GeV^ is 
measured. A cross section value larger than 
the theory predictionJs also reported by the 
ZEUS collaborationEO| in the semi-leptonic 
electron channel. The hh cross section for 
g2 < 1 GeV^ 0.2 < y < 0.8, pt,b > 5 GeV 
and \iib\<2 is: a^p^^,-,^ = 1.6t°:?^nb. The 
NLO prediction is only 0.64 nb. 

A comparison between data and theory 
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Figure 12. Total beauty production cross section in 
77 collisions as a function of the cm energy ^s^e ■ 

for the various measurements is shown in 
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Fig. [llj. For all measurements a higher cross 
section with respect to the NLO predictions 
is measured. However, the total cross sec- 
tion measurement is affected by large extrap- 
olation factors calculated using LO MC sim- 
ulation programs. For instance, jets with 
Et > 6 GeV have to be selected to calculate 
PT,rei ■ In the quoted total cross section they 
are not included. Only 42.9% of the cross 
section is really measured in the detectoica. 
Moreover, requiring two jets with equal Et 
is problematicES. The extianolation includes 
large model uncertaintiesE30 (up to 30%) . It 
is interesting that the ZEUS measurement, 
where also the beauty cross section includ- 
ing the jet requirement is quoted, agrees 
weh with the CASCADE result and with LO 
QCD MC simulations, if heavy quark excita- 
tion processes possibly mimicking part of the 
NNLO correctionscZl are included. 

The increase of the cm energy of LEP 
to y'See ~ 200 GcV has also made possi- 
ble the first observation of beauty produc- 
tion in 7*7* collisions. Here also an excess of 
the measured beauty cross section over the 
theory prediction is founcJ^l (see Fig. |l2| ). 
The NLO prediction agrees with the mea- 
sured charm cross sections, but in the case of 
beauty it underestimates the cross section by 
about two standard deviations. Since many 
of the phenomenological challenges are sim- 
ilar to ep collisions, e.g. the description of 
heavy quarks near the production threshold 
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Figure 13. Beauty production cross section for dif- 
ferent minimal values of the transverse beauty mo- 
mentum in pp collisions at ^/s = 1800 GeV. 

or the treatment of the hadronic structure of 
the photon, this indicates that the excess of 
beauty events at HERA is not accidental, but 
might have a deeper reason. 

Also in pp collision a h igh er beauty cross 
section has been meaauredE3 than expected 
by NLO calculationsEil. In general, the ob- 
served excess grows towards large rapidities, 
i.e. towards the proton remnants. One of 
the measurements performed in the central 
rapidity region is shown in Fig. |l^ over a 
wide minimum transverse momentum range 
p™'". Shown is the b-meson cross section 
(incl. muon) as well as a jet cross section 
where the jets contain a b-tag (muon-|-jet). 
For the last observable the treatment of soft 
and coUinear divergencies is safer in the the- 
oretical calculatioiO. For both observables 
the data lie higher than the theory predic- 
tion. Even so, the agreement with the theory 
is improved in the latter case in particular 
at high p™'"i ^ discrepancy still remains with 
respect to the central result. 

A better agreement with the data is 
found using CASCADE based on the CCFM 
equations^ (see Fig. |l^) . This approach uses 
LO off-shell matrix elements for heavy quarks 
and the scale is set to /^^ ~ ml+p-^ f^. In ad- 
dition, the unintegrated gluon density, which 



has been obtained from an input parameteri- 
sation using only two parameters adjusted to 
HERA data, enters the calculation ^. Note, 
that in a typical NLO DGLAP fit much more 
free parameters are used to parameterise the 
parton densities. In view of the few input 
parameters only fixed by DIS data the agree- 
ment of the CASCADE calculation with the 
data is quite remarkable. A similar success 
in describing beauty production in pp colli- 
sions has been obtained using fcy factorisa- 
tion, an unintegrated parton distribution and 
a BFKL vertex for the gg-vertexEll. 

Whether the higher beauty cross sections 
seen in various processes really can be at- 
tributed to the necessity to modify the par- 
ton evolution at small- or whether it is just 
due to other complications like the resumma- 
tion of large logarithms in the presence of two 
or more hard scales or problems connected to 
the modelling of the fragmentation cannot be 
answered today. The larger amount of data 
which will be available after the HERA and 
TEVATRON luminosity upgrades together 
with the improved efficiency of the upgraded 
detectors will allow further improvements in 
the experimental precision and to extend the 
measurement towards larger pT and larger 
. This might allow possible solutions to 
this problem to be disentangled. 

3.5 Forward Jet and Particle 
Production in DIS 

The classical signature designed to enhance 
BFKL effects Jn DIS is the production of 
"forward jets"E3 characterised by k'^ ~ 
and Xjct = Ejct/Ep as large and x as small 
as kinematically possible (see Fig.||d)|]. The 
first requirement suppresses the kx ordered 
DGLAP evolution. When asking for large 

■'At the TEVATRON cm energy, the gluon-gluon fu- 
sion is the dominant production mechanism. 
* Note that the x values in some of the observables 
are not really small. Moreover, the work to include 
finite 2 terms in CASCADE is still in progress. 
' -Ejct {Ep) denotes the energy of the forward jet (pro- 
ton) in the laboratory system. 
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Figure 14. Forward jet DIS cross section as function 
of X. Approx. NLO BFKL calculations for different 
scales and infra-red cut-offs are shown as lines. 

Xjct/x, the forward jet is separated by a large 
rapidity interval from the struck quark such 
that the phase space for parton emissions be- 
tween the two is amplified. In this case the 
tts In (a;jet/a;) terms are expected to become 
so large that their resummation leads to a siz- 
able increase of the forward jet cross-section: 

ctbfkl oc exp (A In 

X 

A fast rise of the cross-section for forward 
jets with decreasing a; has indeed been seen 
in the dataS. A LO BFKL calculationEl 
rises steeper than the data. A more recent 
calculation emulating some of the expected 
NLO effectsEll is in much better agreement 
(see Fig.[4|). Also CASCADE based lui the 
CCFM equations reproduces the dataO. 

The good description of the forward jet 
cross section can be taken as an indication of 
the need for the BFKL resummation in an ex- 
treme phase space, where small- a; effects are 
expected to be important. However, other 
interpretations are possible. For instance, a 
calculation including direct and resolved pho- 
ton interaction also gives a good description 
of the dataiiil. Kramer and Potter argue 
that the success of their NLO calculation can 
be taken as a strong indication that a direct 
0{al) calculation, i.e. without the need to 
introduce a virtual photon structure, would 
also describe the data. In this case the need 
for a BFKL resummation is just not strong 
enough at HERA to be unambiguously dis- 
tinguished from an exact fixed order (li^alcu- 
lation. However, a recent HI analysis!^ sys- 
tematically investigating the phase space re- 



gions from a simple inclusive jet cross section 
to the forward jet regime shows that neither 
a direct NLO calculation nor the NLO cal- 
culation including resolved photons describes 
the data in the whole phase space. While e.g. 
the forward jet cross section is described, the 
calculation overshoots the simple single inclu- 
sive jet cross section in the backward region. 

In conclusion, despite the wealth of jet 
data at small- a; which became recently avail- 
able a consistent physical interpretation has 
not been achieved. 

A complementary approach to gain in- 
sights in the parton evolution is the produc- 
tion of particles at EtE3. Charged or neu- 
tral particles at high Et are well correlated 
to parton activity. Their production rate 
can be directly predicted by analytic calcula- 
tions folded with proper fragmentation func- 
tions. The measured 7r° cross section exhibits 
the same rise towards small- a; as the inclu- 
sive one, i.e. the ratio Rtt{x) is constant 
(see Fig. At 2 < Q2 < 4 5 Cgy^ 

approximately 0.2% of the DIS events con- 
tain a forward tt" with Et > 2.5 GeV. At 
15 < Q2 < 70GeV^, this ratio rises to 0.5%. 
The MC based on LO QCD and parton show- 
ers (labelled LEPTO) cannot describe this 
behaviour, while a BFKL calculation approx- 
imating some of the NLO effects by requiring 
energy-momentum conservation and allowing 
for a running coupling describes the dataEil. 

4 Inclusive and Diffractitive DIS, 
Transition from Soft to Hard 

The constant (with respect to the inclusive 
cross section) probability to emit hard parti- 
cles in the forward region can be juxtaposed 
to the probability to emit no particle at all. 
In approximately 5 — 10% of all DIS events 
no particle is measured in a large rapidity 
regionED. The latter class of events are called 
"diffractive"0. In both cases the constant 
rate indicates a deep connection between the 
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Figure 15. Rate of forward tt" to the inclusive cross 
section as a function of x. 

production mechanism of rapidity gap events 
and the parton dynamics at small-x. 

The interaction responsible for a large ra- 
pidity gap can be viewed as the exchange of 
a colour singlet decomposing the hadronic fi- 
nal state into a system X and Y: ep eXY 
(see Fig. The kinematic of the process 
can be described by the longitudinal momen- 
tum fraction ^ of the colourless exchange with 
respect to the incoming proton^ and, in anal- 
ogy to Bjorken-x in the inclusive case, by 
the longitudinal momentum fraction (3 that 
is carried by the quark struck by the photon 
(see Fig. |l|). The raticg 



^^(3) 
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is shown in Fig.^ as a function of the 
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In the analogy of the diffraction of light, seq 
refs. therein. 
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"This quanity is often called x^. 
°The diffractive cross section is in this analysis inte- 
grated over My < 1.6 GeV and t < 1 GeV^, where 
t is the squared momentum transfer from the proton 
to the system Y. 




}Mx e = -(l + #) 



9^P /^-I-q™ 



M, 



Figure 16. Sketch of a diflfractive scattering at 
HERA. 



cm energyEa. The ratio p-^^'^^ is relatively 
flat over the full phase space. Only at low 
(3 a tendancy to rise is observed. The con- 
stant rate as a function of the cm energy 
has already been observed in the first HERA 
dataLj, but it is still rather surprising, when 
one tries to understand its physical meaning. 
In the simplest approach^ the colour singlet 
exchange can be modelled by the exchange 
of two "gluons"El. If the "gluons" are inter- 
preted as partons, i.e. pQCD is applicable 
because a hard scale is involved in the pro- 
cess, then at small- a; the inclusive cross sec- 
tion is expected to rise like a oc while 
the diffractive cross section should rise like 
(T cx x^^^ . In this case the rapidity gap 
contribution would significantly increase to- 
wards small-x. If the colour singlet exchange 
is soft, the increase of the diffractive cross sec- 
tion should be closely related to soft hadron 
hadron and 7p collisions where a oc 14^" °^ 
and therefore a oc is expected. In 

this case the rapidity gap contribution in the 
DIS sample would decrease. 

The constant behaviour of p^'-^^ indi- 
cates that the production mechanism of ra- 
pidity gap events contains soft and hard 
pieces. It is remarkable that their interplay 
leads to the same energy dependence of the 
cross section as in the inclusive case. 

For inclusive 7*p collisions, provides 
the hard scale to discriminate soft and hard 
processes. The energy behaviour of the in- 
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Figure 17. Ratio of the diffractive to the inclusive 
DIS cross section as function of the 7*p centre of 
mass energy W in bins of and /3. 

elusive cross section (see eq.^ is shown in 
Fig.^. Recently, ZEUS has achieved a re- 
markable precision on the inclusive cross sec- 
tion measurements at very low thanks 
to dedicated detectors close to the beam 
pipe rneasuring electrons scattered at small 
anglesS For < 0.8 GeV^ a constant 
slope A w 0.1 is found. This value is con- 
sistent with the value (A = 0.08) needed to 
describe the total hadron hadron cross sec- 
tion at high energiesO. The same value also 
describes the energy behaviour of the scat- 
tering of real photons on protonaO. A gen- 
eralised vector dominance modelll3 describes 
the data. For higher values, A linearly 
rises with logQ^ which is well described by 
pQCD. The transition from the soft to the 
hard regime seems to be surprisingly sharp 
(0.5<Q2<1 GeV^). Does this indicate 
that there are really two different physical 
regimes, the long-distance, non-perturbative 
region, where only Regge theory can be ap- 
plied, and the short-distance, pQCD region 
describing the interaction of quarks and glu- 
ons ? Or is a unified description based on 
small-x phenomena possible ? 

The physical understanding of the tran- 
sition region can be made more intuitive in 



the proton rest frame rather than in the DIS 
frame, where the proton has infinite momen- 
tum. In this frame the (virtual) photon splits 
into a qq pair and develops a hadronic struc- 
ture by radiating gluons long before interact- 
ing with the proton. The timescale of this 
fluctuation is r cx l/{xmp). At small-x, the 
lifetime r of the hadronic system is long, e.g. 
at a; = 10~^, cr = 200 fm. The qq pair can be 
viewed as a colour dipole described by a wave 
function 5't,l depending on the longitudinal 
and transverse photon polarisatiorO. It can 
be calculated within QED. The 7*p cross sec- 
tion can be written as: 



cr(a;, 



(frdz |*t,l|^ o-(x,r) 



(6) 



where r is the transverse separation of the qq 
dipole and z is the photon momentum frac- 
tion carried by the quark. By the Heisen- 
berg principle r is closely related to the trans- 
verse quark momentum via r oc l/kr- For 
small kx, i.e. large dipoles, the colour field is 
large and it strongly interacts with the pro- 
ton. The dipole behaves as a hadron in soft 
collisions. For large kx, i.e. small dipoles, 
the colour field is effectively screened and the 
proton is "transparent" to the dipole (colour 
transparency). a{x,r) describes the interac- 
tion of the qq dipole with the proton. 

The following jphenomenological form|^ 
has been proposedL^I for d'{x, r): 



o-Q 1 - exp 



withi?o(x) 



-iIt) 

Xq 



where ctq — 23 mb, xq = 3-10^* and A ~ 0.29 
are determined by a fit to inclusive DIS data 
for X < 0.0l[|. As an example, the dipole 
cross section is shown at = 10 GeV^ in 



Fig. 19, For small r, the dipole cross section 
steeply rises a{x,r) cx r^x~^ as expected by 
pQCD. Towards large r the cross section is 



'Another successful form inspired by the generpljsed 
vector dominance model has been proposed in Ej. 
""The radius Rq has the units 1/GeV. These num- 
bers are quoted for a model with no charm quarks. 
Including charm quarks leads to larger values. 
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Figure 18. Slope of the inclusive total 7p cross sec- 
tion A as a function of Q^. Shown as lines are fits to 
the Golec-Biernat / Wiisthoff model with and without 
QCD evolution. 
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Figure 19. The integrand in equation |6| for inclusive 
(Inc) and diffractive (DD) scattering and the dipole 
cross section a{x,r). 

assumed to saturate a{x,r) oc ctq- It is inter- 
esting that the integrand of equation ^ shown 
in Fig.|l^ peaks around 2/Q making the close 
relation between the photon virtuality and 
the dipole size clear. 

With the assumed form of the fitted 
dipole cross section the inclusive DIS data 
can be successfully described over a wide 
kinematic range (see Fig. |l|b). The addi- 
tional incorporation of DGLAP evolutionll 
is needed to describe the data at large Qy\ 



" To do that, 1/^?q is replaced by 
4:n'^as{fi^)xg{x,fi^)/(3ao) where 



C/r2 + 



Mo- 



However, it can not fully account for the steep 
rise of F2 at large Q^. The transition around 

« 0.5 GeV^ seems to be smoother than 
found in the data. In view of problems tQ_de- 
scribe some of the recent diffractive dataEJ it 
is questionable, if the DGLAP is the correct 
evolution to be applied in this context. 

The diffractive cross section is obtained 
by modifying equation |^ for inclusive DISEZi: 



no 



(frdz |*t,l|^ a'^{x,r) (8) 



By analysing only the leading contribution a 
ratio cx 1 / In {Q^ Rl{x)) is found 

which is only slowly varying with x and 
and correctly reproduces the data0 when 
in the dipole also qqg states are included. 
Fig. ^ shows that the integrand of eq. || 
contains a significantly higher contribution of 
large dipole sizes at the same than found 
in the inclusive case. This makes both the 
soft and the hard nature of the diffractive 
production mechanism evident. By assum- 
ing a certain interaction of the colour dipole 
with the proton based on the inclusive DIS 
data at small-x the Golec-Biernat/Wiisthoff 
model predicts the diffractive cross section. 
This underlines the deep connection between 
diffraction and the small- a; parton dynamics 
which we start to understand better and bet- 
ter. In the leading logarithm approximation 
the colour dipole formulation and the kr fac- 
torisation give an equivalent description of 
hard processes at high energjnia. 

5 Instanton Production in DIS 

In QCD, anomalous non-perturbative pro- 
cesses are expected which violate classical 
laws such as Jjae conservation of chiral- 
ity. InstantonsEJ, non-perturbative fluctua- 
tions of the gluon field, represent tunnelling 
transitions between topologically inequiva- 
lent vacua. DIS offers a unique possibility 
to discover QCD instantons. 




Figure 20. Sketch of the production of instanton in- 
duced events in DIS. 

In DIS instanton induced processes (J) 
are dominantly produce d iu a, quark gluon fu- 
sion process (see Fig.|^) E§Ea. The virtuahty 
Q'^ of the quark g', originating from a photon 
splitting into a qq pair in the /-background, 
provides a generic hard scale naturally lim- 
iting the /-size p and makes a quantitative 
prediction of the cross-section possible eSeZI. 
The expected cross section for the kinemat- 
ical region x > 10^^ and (LI <; y < 0.9 is 
of the order of 10 - 100 pbia@. This cross 
section is sizeable, but still much lower than 
the inclusive DIS cross section. The expected 
characteristics of the hadronic final state of /- 
processes, which can be simulated using the 
QCDINSE3 program, can be exploited to im- 
prove the signal to noise ratioEll. Besides the 
current jet, a large number of particles of all 
kinematically allowed flavours is produced at 
high Et in a narrow band of two rapidity 
units. Since the Instanton decays isotropi- 
cally, a large sphericity of the particles in the 
/-rest system is expected. 

The _H1 collaboration has recently 
performedES the first dedicated search for 
QCD /-induced processes using six hadronic 
final state /-sensitive observables in the kine- 
matic range x > 10~^, 0.1 < y < 0.6 and 
6ei > 156°, where 6^1 is the polar angle of 
the scattered electron^ As an example the 
sphericity distribution is shown in Fig. [2l| a 
for the inclusive DIS sample. Before cuts to 
enhance the /-fraction, the data are well de- 
scribed by pQCD models. The expected in- 



xg{x,fj.^) = AgX ^9 and C, Ag, Xg and fi'^ are addi- 
tional fit parameters. 



'For an /-analysis based on existing data see | 
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Figure 21. Sphericity for normal and /-DIS without 
(a) and with (b) cuts to enhance the /-signal. 

Stanton sample is about 2-3 orders of mag- 
nitude smaller and the events are expected 
to be more spherical. After cuts the back- 
ground is suppressed by about a factor of 
1000. 549 events are selected in the data, 
while only 363126 ^^'^ 435^22 events are ex- 
pected by two different pQCD models. The 
size of the excess is, however, about the same 
as the difference in the pQCD background 
models. The sphericity distribution of these 
events is shown in Fig.^l|b. In four out of 
the six /-sensitive observables a qualitative 
similarity between the difference of the data 
and the pQCD models and the expected I- 
shape is observed. The two remaining ob- 
servables do not particularly support the hy- 
pothesis that the observed excess can be ex- 
plained by I-processes as currently imple- 
mented in QCDINS. However, recently it has 
been pointed out that the shape of these two 
particular observables is affected by large the- 
oretical uncertainties due_to missing cuts in 
the experimental analysia^J. The data limit 
the allowed /-cross section to about 100 pb. 
A steep rise of the I-cross section towards 
large instanton sizes, as would be obtained 
from a naive extrapolation of /-perturbation 
theory, is excluded by the data. The absence 
of such a steep rise is in agreement with lat- 
tice simulations of zero flavourfl'S. 

Conclusions 

The strong cross section increase expected in 
LO BFKL for various processes has not been 
observed. Thanks to a huge theoretical effort 
significant progress has been recently made. 
The BFKL NLO corrections have been calcu- 
lated and additional resummations have been 



worked out. The results are in qualitative 
agreement with the data. However, no clear 
and unmistakable evidence for the need for 
small-x effects beyond the standard resum- 
mation of leading logarithms using coUinear 
factorisation has been reported so far. Nev- 
ertheless there are strong indications that 
small-x effects play already at the present en- 
ergies a certain role and that they describe 
high energy collisions in a more coherent way. 
First encouraging results using angular order- 
ing as a key unifying principle have been ob- 
tained. This allows more precise QCD pre- 
dictions for future colliders to be made. 

The increased precision of the HERA 
measurements reveals more and more the in- 
triguing connection between soft and hard 
physics. In the understanding of soft phe- 
nomena the correct description of diffraction 
plays a key role. Here, the deep connection 
to small-x parton dynamics becomes clearer 
and clearer and more is to be learned in the 
near future. The final aim is to get bet- 
ter insights to the confinement of quark and 
gluons in hadrons. The pioneering experi- 
mental and phenomenolgical work to discover 
QCD instanton-induced processes might add 
important complementary information in our 
understanding of non-perturbative phenom- 
ena. It seems possible to establish such an 
important non-perturbative effect at HERA, 
although no clear experimental evidence has 
been found so far. 
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Questions: 

G. Altarelli, CERN: 

The beautiful HERA data on small- a; struc- 
ture functions pose a clear problem: Why do 
the NLO QCD evolution equations work so 
well in spite of the fact that the BFKL correc- 
tions could be large. You have mentioned the 
Salam, Ciafaloni, Colferai proposed solution 
to this problem, but nobody has so far shown 
that their theory fits the data well. Myself, 
Ball and Forte we have developed an alter- 
native approach, less predictive, but more 
model independent (in our opinion) which we 
have proven to fit the data (we can accom- 
modate but not predict the smallness of the 
BFKL corrections). Until the understanding 
of totally inclusive structure functions is not 
more satisfactory, the theoretical foundations 
of many BFKL signals in non inclusive distri- 
butions is quite shaky. 
Answer: 

I agree that the F2 data can be interpreted 
in the theoretically cleanest way, since the 
uncertainties are usually much smaller than 
for less inclusive observables. I have pointed 
out that many open questions remain in the 
understanding of F2 and I have stressed the 
importance of even more precise data. Since 
F2 is a very inclusive quantity, precision is the 
key issue. However, less inclusive observables 
allow to look at new effects in a more direct 
way and therefore might be more promising. 
At the end of the day, it will be important 
to have a consistent theory of strong interac- 
tions to describe inclusive and less inclusive 
measurements in many different processes. 

P. Nason, INFN Milano: 

I have a remark with regard to the CCFM 
Monte Carlo predictions. Monte Carlo pro- 
grams are practical implementation of theo- 
retical ideas and calculations, and at the end 
they should match precise analytical calcu- 
lations. For example, attempts to explain 
the beauty cross section by small- a; resum- 
mations, once 0{a1) and 0{a^g) corrections 



are correctly implemented, failed (until now) 
to explain the data. Thus, the inaccuracy 
of the Monte Carlo (that certainly does not 
implement the ©(a^) correctly) should be 
addressed before claims are made of solving 
phenomenological problems. 
Answer: 

The CASCADE MC is an implementation 
of the CCFM equation. It seems to be ad- 
equate, since in all proposed CCFM solu- 
tions the correct treatment of the kinematics 
turned out to be important. Although it is 
based on LO, part of the higher orders are 
resummed. In particular when factorisa- 
tion is used, it is difficult to say if, e.g. a 
©(a^) calculation based on the coUinear fac- 
torisation contains more or less higher order 
contributions. In this sense there is no inac- 
curacy, but it is a different approximation. 

K. Ellis, Fermilab: 

You showed calculations due to Jung and 
Salam which lead to a larger beauty quark 
cross section at the TEVATRON. It is im- 
portant to establish what feature of the cal- 
culations give rise to the increase. Does the 
unintegrated gluon distribution used in the 
calculations give (after integration over fcf.) 
a gluon distribution in accord with the gluon 
distribution measured elsewhere at HERA ? 
Answer: 

The gluon density is not a physical ob- 
servable. Therefore the gluon as extracted 
from NLO DGLAP fits and from the CCFM 
scheme do not have to agree and in fact they 
do not. The unintegrated gluon is higher Eil. 
However, it is only important that both de- 
scribe the cross section measurements and in 
fact they both do so. 
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